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Abstract 

We present dynamical description of gravitational collapse in view 
of Misner and Sharp's formalism. Matter under consideration is a 
complicated fluid consistent with plane symmetry which we assume 
to undergo dissipation in the form of heat flow, radiation, shear and 
bulk viscosity. Junction conditions are studied for a general spacetime 
in the interior and Vaidya spacetime in the exterior regions. Dynam- 
ical equations are obtained and coupled with causal transport equa- 
tions derived in context of Miiller Israel Stewart theory. The role of 
dissipative quantities over collapse is investigated. 
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1 Introduction 

The ultimate fate of the star (when it undergoes catastrophic phase of col- 
lapse) is one of the most important questions in gravitation theory today. 
When a star has exhausted all of its nuclear fuel, it collapses under the in- 
fluence of its own gravity and releases large amount of energy. In fact, it is 
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a highly dissipative process, i.e., energy is not conserved in it, rather due to 
various forces and with the passage of time, it becomes lesser. Dissipative 
process plays dominant role in the formation and evolution of stars. 

The initial discussion over this problem was given by Oppenheimer and 
Snyder [1] who assumed a spherically symmetric distribution of matter. They 
took the most simplest form of matter, i.e., dust and the flow is considered 
to be adiabatic. It is somewhat unrealistic to ignore the pressure as it cannot 
be overlooked in the formation of singularity. Misner and Sharp adopted 
a better approach by considering an ideal fluid which gave a more realis- 
tic analysis of gravitational collapse. Both of them assumed vacuum in the 
exterior region. Vaidya [3] introduced a non-vacuum exterior by giving the 
idea of outgoing radiation in collapse. It was physically a quite reasonable 
assumption as radiation is a confirmation that dissipative processes are oc- 
curring, causing loss of thermal energy of the system which is an effective 
way of decreasing internal pressure. 

The Darmois junction conditions [1] gave a way to obtain exact models of 
an interior spacetime with heat flux to match with exterior Vaidya spacetime. 
Sharif and Ahmad [S] considered the perfect fluid with positive cosmological 
constant to discuss the junction conditions with spherical symmetry. The 
same authors [6] also worked on junction conditions for plane symmetric 
spacetimes. 

Goswami [7] made an attempt in search of a more physical model of 
collapse. He considered dust like matter with heat flux to conclude that 
dissipation causes a bounce in collapse before the formation of singularity. 
Nath at el. [8] investigated dissipation in the form of heat flow and formu- 
lated junction conditions between charged Vaidya spacetime in exterior and 
quasi-spherical Szekeres spacetime in interior regions. They also discussed 
apparent horizons and singularity formation. Ghosh and Deshkar [9j studied 
gravitational collapse of radiating star with plane symmetry and pointed out 
some useful results. A lot of work is being done over gravitational collapse by 
considering shear free motion of the fluid. Although, it leads to simplifica- 
tion in obtaining exact solutions of the field equations, yet it is an unrealistic 
approach. Shear viscosity is a source of dissipating energy and plays an im- 
portant role in collapse. Chan [10] investigated gravitational collapse, with 
radial heat flow, radiation and shear viscosity. He showed how the pressure 
became anisotropic due to shear viscosity. 

Herrera and Santos [llj discussed the dynamics of gravitational collapse 
which undergoes dissipation in the form of heat flow and radiation. Di Prisco 
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et al. [12] extended this work by adding charge and dissipation in the form 
of shear viscosity. Herrera [IB] provided comprehensive details of inertia of 
heat and how it plays an effective role in dynamics of dissipative collapse. 
Herrera and Martinez [H] presented relativistic model of heat conducting 
collapsing object and debated over the effect of a parameter which occurs 
in dynamical equation on collapse. Herrera and collaborators pro- 
posed a model of shear free conformally flat collapse and focused on the role 
of relaxation process, local anisotropy and relation between dissipation and 
density inhomogeneity. 

Recently, Herrera et al. [17] threw light on behavior of non-equilibrium 
massive object which lost energy due to heat flow, radiation, shear and bulk 
viscosity. Matter under consideration was distributed with spherical symme- 
try. It has become quite clear that when mass and energy densities involved 
in the physical phenomenon are sufficiently high as in gravitational collapse, 
gravitational field plays an important and dominant role. The gravitational 
dynamics then must be taken into account for a meaningful description of 
such ultra high energy objects. This fact motivated us to elaborate the 
above mentioned paper in the context of plane symmetries. Matter under 
consideration is a complicated fluid which suffers through dissipation. Misner 
and Sharp's prescription is used to work out dynamical equations. Transport 
equations are obtained in the context of Miiller Israel Stewart theory [18], [19] 
which is a causal theory for dissipative fluids. Thermodynamic viscous/heat 
coupling coefficients are taken to be non-vanishing which is expected to be 
quite plausible in non-uniform stellar models of universe. One of the dynam- 
ical equations is then coupled to transport equations in order to figure out 
the influence of dissipation over collapse. 

The paper is written in the following manner. The next section is about 
the matter distribution in the interior region and some physical quantities rel- 
evant to matter under consideration. The Einstein field equation are worked 
out in section 3 and junction conditions are discussed in section 4. Dy- 
namical equations are formulated in section 5 and are coupled to transport 
equations in section 6. The last section discusses and concludes the main 
results of the paper. 
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2 Interior matter distribution and some phys- 
ical quantities 

A 4-dimensional spacetime is split into two regions: interior y~ and exterior 
through a hypersurface S which is the boundary of both regions. We 
assume the matter distribution in the interior region to be consistent with 
plane symmetry. The interior region V~ admits the following line element 

ds'i = -f{t, z)dt^ + git, z){dx^ + dy"^) + h{t, z)dz^, (1) 

where {x~'^} = {t, y, z} (/i = 0, 1, 2, 3). The fluid is presumed to dissipate 
energy in terms of heat flow, radiation, shearing and bulk viscosity. 
The energy-momentum tensor for such a fluid is defined as 

Tab={fl + P + n)KH +{P + ^)9ab + QaVb + ^feK + djb + T^ab, (2) 

where /i, p, U, qa, la and are the energy density, pressure, bulk viscosity, 
heat flow, null four-vector in z-direction and shear viscosity tensor respec- 
tively. Heat flow Qa is taken to be orthogonal to velocity V"", i.e., QaV"" = 0. 
Moreover, we have 

^"¥^ = -1, rva = -i, 71,^^ = 0, iT[ab] = o, < = o, n, = o. (3) 

In the standard irreversible thermodynamics by Eckart, we have the following 
relation [20] 

TVab = -2r]aab, n = -Ce, (4) 

where r] and C stand for coefficients of shear and bulk viscosity, aab is the 
shear tensor and G is the expansion. The algebraic nature of Eckart constitu- 
tive equations causes several problems but we are concerned with the causal 
approach of dissipative variables. Thus we would not assume (jl]) rather we 
shall resort to transport equations of Miiller-Israel- Stewart theory. 
The shear tensor aab is defined as 

(Tab = V(a;6) + a(aH) " ^©/iafe, (5) 

where the acceleration aa and the expansion B are given by 

aa = Va;bV\ e = V;, (6) 
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and hab = gab + Vayb is the projection tensor. The shear tensor Gab satisfies 

V;a<^^' = 0, a^^ = a^^, < = 0. (7) 
In co-moving coordinates, one can take 

^°^> '°^77'» + 7s*- 

here g is a function of t and z. 

Using Eq.([8]), the non- vanishing components of the shear tensor aab turn 
out to be 

9 2/i 
= --0- = 0-22, cr33 = — cr, (9j 
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where 



Thus we have 



,= J_f^_9V (10) 

9 

O^afcO- = -cr . (llj 

Also, in view of Eqs.(l3]) and (jl]), it yields 

7roa = 0,7r| = -27r2 = -27rJ. (12) 
In compact form, it can be written as 

TTafe = fi(XaXb - ^^afe), (13) 

where f2 = and x*^ is a unit four-vector in z-direction satisfying 

X"Xa = l, X"K = 0, X'' = ^^3- (14) 
In view of Eqs.(l6]) and (jH]), it follows that 

«3 = 4, = ^f^ + 4:V (15) 

where dot and prime represent derivative with respect to time t and z re- 
spectively. 

The Taub's mass for plane symmetric spacetime is defined by [21] 
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3 The Einstein field equations 

The Einstein field equations for tlie metric ([1]) yield the following set of 
equations 

9 I 9 , h\ fg' (h' g'\ fg" 



2g \ 2g h I 2gh \h 2g J gh 



2f \2f 2g 2hlAh\f h g J Afh\h f 

(18) 

9 /«' , /'^ , 9'h 9' + (20) 



2fi' Vfi' // 9 
After some manipulation, we can also write Eq.(l20l) in the following form 

47r(g + e)Vh = \{Q-a)' - a^. (21) 

4 Junction conditions 

We discuss junction conditions for the interior region V' given by Eq.([T]) and 
the exterior region which is taken as plane symmetric Vaidya spacetime 
ansatz given by the line element |22] 

dsl = ^^du'' - 2dpdZ + Z\dX^ + dY^), (22) 

where = {^^ ^^Y, Z} (/i = 0,1,2,3), v is the retarded time and m{v) 
represents total mass inside S. The line element for the hypersurface S is 
defined as 

[ds^)^ = -dr^ + A^{T){dx^ + dy"^), (23) 
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where ^* = (r, x, y) {i = 0, 1, 2) are the intrinsic coordinates of S. 
The Darmois junction conditions |1] are 

• The continuity of the hue elements over the hypersurface S gives 

{ds')j: = (rf^)E = {dslh. (24) 
This is called continuity of the first fundamental form. 

• The continuity of the extrinsic curvature Kab over the hypersurface E 
yields 

[Kij]=Kt,-Kr: = 0, (a, 6 = 0,1, 2). (25) 
This is known as continuity of the second fundamental form. 

Here K^i is the extrinsic curvature defined as 



fS = -"J(a^ + r;;„^L_^). (.,^. = 0,1.2.3). (26) 



where are the components of outward unit normal to hypersurface E in 
the coordinates x^^- 

The equations of hypersurface E in terms of coordinates x^'^ are given as 

k^t,^) = z-zj: = 0, (27) 
k4u,Z) = Z-Zj,{u) = 0, (28) 

where z^, is taken to be an arbitrary constant. Using Eqs. fl27|) and fl28|) . the 
interior and exterior metrics take the following form over hypersurface E 

{dsl)Y: = -f{t,zj,)dt^ + g{t,zj:){dx^ + dy^), (29) 
2^ n( dZj: 



(dsih = 2\^-^~^jdu' + ZUdX' + dY'). (30) 

In view of junction condition ( l24l) . we get 

Zl = g{t,zj,), (31) 
dt 1 



dr v7' 

di^ 1^2^^'^ 2m (z/) 

dr \ du Zy. 
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-1/2 



(32) 
(33) 



Using Eqs.f l27p and f l28l) . the unit normals in V and respectively, turn 
out to be 



n, 



n 



v^(0, 0,0,1), 
\du Z 



-1/2 



dZ 



(34) 
(35) 



The non-zero components of the extrinsic curvature are 



00 



K, 



00 



d^u f dv\ ^ m dv 
d^\Jh) ~~Z^d^ 



J s 



K 
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K. 



22 



^^2^2 



9' 



2Vh 



Zf- - 2m^ 
dr dr 



(36) 

(37) 

(38) 
(39) 



Now, by the junction condition f l25|) . i.e., continuity of extrinsic curva- 
tures, it follows that 



d'^u / di'\ ^ m dv 
d^Uh) "Z^d^ 



2m 



dp g 



9' 



dr 2v7 2Vh 
Using Eqs.( l33|) and (14T]) . we obtain 

f dv\ ^ 



9^9 



2v7 2Vh. 



(40) 
(41) 

(42) 



Inserting Eq.(j42]) in dH]), it follows that 



m[v) 



1 fg'^ 9' 



8v^ V/ h 



and hence 



m(t, z) = m(u). 



(43) 
(44) 
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Differentiating Eq. (H^ with respect to r, and making use of Eqs. (H5il and 
we can write Eg. (HQ]) as 



"/h j -9 gf 



-g' ^ g'h ^ f'g ^ 



f (^\ ' ( g'g 



0. (45) 



Comparing Eq.(j45]) with Eqs.dlH]) and ([20]), it yields 



p + n + -n = g. 



(46) 



5 Dynamical equations 

The energy-momentum conservation, = 0, gives 



+ 



(ii + e) (q' + e') g , ^ 1^, 



and 



e + fi) = 0. (48) 



Now we investigate the dynamical properties of the system using the Mis- 
ner and Sharp's [2] perspective. For this purpose, we take the proper time 
derivative as 

and the proper derivative in z-direction as 

1 d 



Z'dz' 



(50) 



where 



Z = ^. (51) 
The velocity U of the collapsing fluid can be defined as the variation of 



Z with respect to the proper time 

U = Dt{Z) 



1 



DtQ. 



(52) 



In the case of collapse, the velocity of the collapsing fluid must be negative. 
In view of Eq.(l52l). Eq.( lT6l) can take the following form 



E 



Making use of Eq.dSO]) in Eq.dH]), it follows that 
A-K^q + e) = E 



5^.(8-.)-! 



In case of no dissipation, using Eqs. lfTO]) . ( IT^ and 
becomes 

= 0. 



(53) 
(54) 

the above equation 
(55) 



This implies that f/ ~ Z depicting that now collapse will be homologous. 
The rate of change of Taub's mass, using Eqs. f|T6|) . f|T9|) . fl20|) and fH9|) . turn 
out to be 

DtTU = -AnZ^lip + n + e + -n)U + {q + e)E]. (56) 

o 

Thus the rate of change of Taub's mass represents variation of total energy 
inside the collapsing plane surface. Since this variation is negative, it shows 
that total energy is being dissipated during collapse. The first round brackets 
on the right hand side stand for energy due to work being done by the effective 
isotropic pressure {p + U + |r2) and the radiation pressure e. The second 
brackets describe energy leaving the system due to heat flux and radiation. 
Similarly, using Eqs.([IS]), ([IT]), 020]) and (EU]), we get 



D^m = 47rZ^[// + e + {q + e 



E^ 



(57) 



This equation describes about the variation of energy between adjoining plane 
surfaces inside the fluid distribution. On the right hand side, (// + e) stands 
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for energy density of the fluid element plus the energy of null fluid showing 
dissipation due to radiation. Moreover, {q + e)^ is negative (as U < 0), 
telling that energy is leaving due to outflow of heat and radiation. 

Making use of Eqs.(Il6]), (Il9]), (EI]) and ([53]), the acceleration DtU of the 
collapsing matter inside the hypersurface S is given as 

3 ^ Z2 2fVh ^ ' 

Substituting the value of ^ from the above equation into Eq. fHS]) . it follows 

that 

(p + n + /i + 2e+ ^^])Dr^ = -(p + n + /x + 2e+ ^^]) 

x[47rZ(p + n + e + -fi) + ^] 

-E\D^[p + n + e + ^^]) + |(e + ^)\ 

-EIDtQ + Dre + 4(g + e)£ + 2{q + e)a]. 

(59) 

This equation has the form of Newton's second law, i.e.. 

Force = Mass density x Acceleration. 

The term within the brackets on the left hand side stands for "effective" 
inertial mass and the remaining term is acceleration. The first term on 
the right hand side represents gravitational force. Since by the equivalence 
principle, inertial mass is equivalent to passive gravitational mass and passive 
gravitational mass is equivalent to active gravitational mass. Thus the factor 
within round brackets stands for active gravitational mass and the factor 
within the square brackets shows how dissipation effects active gravitational 
mass. The second square brackets firstly include gradient of effective pressure 
which involves radiation pressure and the collective effect of shear and bulk 
viscosity. The second contribution is of local anisotropy of pressure which is 
the result of radiation and shear viscosity. The last square brackets entirely 
depend upon dissipation. The hydrostatic equilibrium can be obtained from 
the above equation by substituting U = 0, q = 0, e = 0, 11 = and i7 = 0. 



DzP= -ifJ' + P)-^ 
Z' 



ITT- , ^ 

+ A-kZp 

Z2 
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6 Transport equations 



The general expression for entropy 4-current is given as [20 



2T 



+ 



T 



+ ■ 



T 



-, (60) 



where n is particle number density, T is temperature, /3yi(p, n) > are ther- 
modynamic coefficients for scalar, vector and tensor dissipative contributions 
to the entropy density and ayi(p, n) are thermodynamic viscous/heat cou- 
pling coefficients. The divergence of extended current (follows from Gibbs 
equation and Bianchi identities) is given by 



-n 



n 



- g"K(lnr),^(l + aoH) + V^,^V^ - a^U,^ - ai<.^ 



IT 



cr, 



T f P2 



(61) 



The 2nd law of thermodynamics requires that S"." > 0. This leads to the 
following transport equations for our dissipative variables 



ron„y" + n = -ce + aoCCa-^CT(^^v^°) n. 



(62) 



rihiqp.^V^' + q^ = -A;[/i^T^(l + aoH) + ai^/ij^T^ + T(a„ 

- aoH;. - ai<^)] - hr' {^V^) g« (63) 



and 



T2h'^hp7r^„.pV^ + tTq/s = -2rjaap + 2r]aiq^p.a> - vT TTTf^" '^ap, (64) 



T2 



where 



2r]T 



(65) 
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with k as the thermal conductivity. The relaxation times are given by 



To = C^o, n = kTjSi, T2 = 2rij32. 



(66) 



Notice that if the thermodynamic coupling coefficients are assumed to be 
zero, Eqs.(|62D-([64D turn to be Eqs.(2.21)-(2.23) as given in [20]. The inde- 
pendent components of Eqs.( l62|) -( |64|) are calculated as follows. 



2 \CT 



n. 



(67) 



-k 



v7 
2 



T'(l + aoU + ^al^]) + T - aoll" 
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(68) 



Vh g 



r]T 



r2 

2r]T 



(69) 



Now we discuss the action of dissipation over dynamics of collapsing ob- 
ject. We couple these transport equations to dynamical equation (j59|) . Using 
Eq.(l6HD in Eq.(l59]), it follows that 



(/i + p + n + 2e + -(])(! - A)DtU = (1 - A)Fgrav + F^yd 



kE 



2 r 



kT'^q 
~2^ 



3 



- E 



3q 



u 



^ + 2e e - - - 2(g + e)^ 



Z 



(70) 



13 



where Fgrav and Fhyd are given by 



and 



grav 



-{p + U + fi + 2e + -n) 

o 

2 

m + 47r(p + n + e + -Q)Z'^ 
3 



1 

i2' 



hyd 



D^(p + n + e + ^fi)+2(e + l])i 



kT 



A = — p + n + ^ + 2e+-fi 



1 - -aiQ 



Inserting Eg. (1671) in Eq. (170|) . we obtain 



(71) 
(72) 

(73) 



(p + n + /i + 2e + -^])(l - A + A)DtU = (1 - a + A)Fgrav + Fhyd 
kE"^ r / 2 \ ( 2 / 3 

^ + 2(g + e)£' 



-E' [p + n + ti + 2e + -n] A 



IkT^q f Tl \ 
A 



+ E 



+E 



p + ^ + /i + 2e + -^] 



n + toDtU 



(74) 



where A is given by 



aoCq {p + U + n + 2e + -n 



3g + 4e 
2C + Ton 



(75) 



Here we see that (1— A+A) is the major factor that appears in the dynamical 
equation after coupling it with the transport equations. We would like to 
mention here that Eq. flT^ is the plane symmetric version of Eq.(55) in [IT]. 



7 Summary and Conclusion 

Gravitational collapse in a star is an irreversible phenomenon. Dynamics 
(such as transport processes) of such non-equilibrium objects and connection 
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between their dynamics and thermodynamics are of extensive significance in 
order to have a better visuahzation of this problem. Thus we have studied 
the dynamics of dissipative collapse, i.e., what role does dissipation play 
with passing time as star collapses under the infiuence of its own gravity. 
The most realistic model of matter, i.e, complicated fiuid is assumed in the 
interior region and is taken to be consistent with plane symmetry. 

To see how system evolves with time, dynamical equations for the plane 
symmetric spacetime are obtained using Misner and Sharp formalism. In 
the dynamical equation (!59|) . we see that the gravitational force represented 
by the first term on the right hand side is expected to be much effective as 
compared to non-dissipative fiuid and so gravitational collapse is expected to 
be faster in this case. Moreover, since the pressure gradient is negative in the 
second term on right hand side of this equation, which combined with the mi- 
nus sign preceding that term makes a positive contribution, thereby reducing 
the rate of collapse. The last square brackets entirely depend on dissipation 
and one cannot expect any such contribution in a dynamical equation for 
non-dissipative collapse. The third term in this bracket is positive due to 
negative sign of velocity of collapsing fiuid U. It shows that outfiow of heat 
fiux g > and radiation e > reduces the total energy of the system and 
hence reduces the rate of collapse. 

Transport equations in the context of Miiller, Israel and Stewart theory 
of dissipative fiuids are obtained and coupled to dynamical equation in order 
to see the infiuence of dissipation over dynamics of a collapsing plane. After 
this union of dynamical and transport equations, we get equation fl7^ where 
the factor (1 — A + A) appears in the dynamical equation. We see the effect 
of this factor for different possible values. 

• IfO<(A — A)<1, inertial and gravitational mass densities will be 
reduced. 

• If (A — A) tends to 1, inertial mass density tends to zero. 

• If (A — A) > 1, gravitational force will become positive and it will lead 
to the reversal of collapse. Another possibility for reversal of collapse 
is to take (A — A) < 1 such that (1 — A + A) is sufficiently small. 
Consequently, it will significantly decrease the gravitational force. 
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